Zeno and anti-Zeno effects for quantum Brownian motion 
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In this paper we investigate the occurrence of the Zeno and anti-Zeno effects for quantum Brownian 
motion. We single out the parameters of both the system and the reservoir governing the crossover 
between Zeno and anti-Zeno dynamics. We demonstrate that, for high reservoir temperatures, 
the short time behaviour of environment induced decoherence is the ultimate responsible for the 
occurrence of either the Zeno or the anti-Zeno effect. Finally we suggest a way to manipulate the 
decay rate of the system and to observe a controlled continuous passage from decay suppression to 
decay acceleration using engineered reservoirs in the trapped ion context . 
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The quantum Zeno effect (QZE) predicts that the de- 
cay of an unstable system can be slowed down by measur- 
ing the system frequently enough Q. In some systems, 
however, an enhancement of the decay due to frequent 
measurements, namely the anti-Zeno or inverse Zeno ef- 
fect (AZE), may occur 0- 

In this paper we focus on the quantum Brownian mo- 
tion (QBM) model 0,110 which is a paradigmatic model 
of the theory of open quantum systems. This model, 
dealing with the linear interaction of a particle with a 
bosonic reservoir in thermal equilibrium, is widely used 
in several physical contexts. It describes the dynamics of 
a particle interacting with a quantum field in the dipolc 
approximation as well as a quantum electromagnetic 
field propagating in a linear dielectric media 7]. The 
model is used in nuclear physics to describe, e.g., the 
two-body decay of an unstable particle ||. In quan- 
tum chemistry, the QBM model describes the quantum 
Kramers turnover, which forms the basis of modern the- 
ory of activated processes || . 

Recently, this model has been investigated to explain 
the loss of quantum coherence (decoherence) due to the 
interaction between the system and its surroundings. In 
particular, the absence of macroscopic quantum super- 
positions in the classical world has been explained, using 
the QBM model, in terms of environment induced de- 
coherence (EID) jfjj. With the last term we mean here 
a process which transforms a highly delocalized state in 
position and/or momentum, e.g. a superposition of co- 
herent states, into a localized classical state. 

To the best of our knowledge, the conditions for the oc- 
currence of the quantum Zeno and anti-Zeno effects have 
never been considered for the QBM model. The aim of 
this paper is to investigate the Zeno and anti-Zeno phe- 
nomena in a system, namely the damped harmonic os- 
cillator, which possesses a classical limit and where it is 
therefore possible to monitor the transition from quan- 
tum to classical dynamics caused by decoherence induced 
by the environment. Previous studies focus on few level 
systems and deal with completely quantum states, such 
as the spin, which have no classical analogue [lfj, . 

Our main result is the demonstration that the occur- 
rence of either the Zeno or the anti-Zeno effect stems 



from the short time behaviour of the environment in- 
duced decoherence, which therefore drives the Zeno-anti- 
Zeno crossover. On the other hand, one can use the QZE 
or the AZE to manipulate the quantum-classical border 
by prolonging or shortening, respectively, the persistence 
of quantum features in the initial state of the system. 
Moreover, we suggest a physical context in which the 
Zeno-anti-Zeno crossover can be observed with current 
technology by means of reservoir engineering techniques 

The QBM microscopic Hamiltonian model consists of a 
quantum harmonic oscillator linearly coupled to a quan- 
tum reservoir modelled as a collection of non-interacting 
harmonic oscillators at thermal equilibrium. In the limit 
of a continuum of frequencies u>, the reservoir proper- 
ties are described by the reservoir spectral density J(oS) 
measuring the microscopic effective coupling strength be- 
tween the system oscillator and the oscillators of the 
reservoir. 

One of the advantages of the QBM model is that, for 
factorized initial conditions, it can be described by means 
of the following exact master equation 0, 0, 1 1 ll ] 
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where ps(t) is the reduced density matrix, and Hq = 
fkvo (a^a + |) is the system Hamiltonian, with a, a', 
and ujq the annihilation operator, the creation oper- 
ator and the frequency of the system oscillator, re- 
spectively. Moreover X = yUOo/2h (a + at) and P = 
i/-\/2Tujj( S [a/ — a) are the system position and momen- 
tum operators, where we have set the mass of the particle 
m = 1. 

The master equation given by Eq. Q), being exact, 
describes also the non-Markovian short time system- 
reservoir correlations due to the finite correlation time 
of the reservoir. In contrast to other non-Markovian dy- 
namical systems, this master equation is local in time, 
i.e. it does not contain memory integrals. All the non- 
Markovian character of the system is contained in the 
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time dependent coefficients appearing in the master equa- 
tion (for the analytic expression of the coefficients see, 
e.g., Ref. Q). These coefficients depend uniquely on the 
form of the reservoir spectral density. The coefficient 
r(t) describes a time dependent frequency shift, j(t) is 
the damping coefficient, A(t) and H(t) are the normal 
and the anomalous diffusion coefficients, respectively 
In the secular approximation, i.e. after averaging over the 
rapidly oscillating terms appearing in the dynamics, the 
only two relevant coefficients are the diffusion coefficient 
A(t) and the damping coefficient 7(f) A simple view 
of the effect of the diffusion term can be given by looking 
at the approximate solution Q of the master equation 
given by Eq. 0J, in the position space, 



p s (x, x', t)~p s {x, x', 0) exp 



-(x- 



x') 2 Ait^dh 
Jo 



(2) 



The previous equation shows that the diffusion term is 
responsible for the vanishing of the off-diagonal terms of 
the density matrix in position space, i.e. for EID. 

For the sake of concreteness, we focus on the case of 
an Ohmic spectral density with Lorentz-Drude cutoff (see 

0, P-25), J(w) — (llj/tt)u! 2 /(llj 2 +u! 2 ) where lu c is the cut- 
off frequency. This form of the spectral density is one of 
the most commonly used since it leads to a friction force 
proportional to velocity, which is typical of dissipative 
systems in several physical contexts. The main result of 
the paper, however, holds for general forms of the spec- 
tral density. 

We assume that the system oscillator is initially pre- 
pared in one of the eigenstates of its Hamiltonian, i.e. a 
Fock state |n). During the time evolution the system 
is subjected to a series of non-selective measurements, 

1. e. measurements which do not select the different out- 
comes 0| • We indicate with r the time interval between 
two successive measurements, and we assume that r is 
so short (and/or the coupling so weak) that second order 
processes may be neglected. Stated another way we as- 
sume that the probability P n {t) of finding the system in 
its initial state \n), i.e. the survival probability, is such 
that P„(i) ~ 1 and P n (t) > P„±i(t). After TV measure- 
ments the survival probability reads [ToL Hr3 | 



pW(t)=P n (r) If = eX p [-^{r)t] 



(3) 



where t = Nt is the effective duration of the experiment 
and the effective decay rate 7„ (r) is defined by the last 
equality. In Eq. ©, we have assumed that the proba- 
bility P n (r) factorizes. This assumption is justified by 
the fact that, in the following we will use second order 
perturbation theory. As shown, e.g., in Ref. 17], up to 
second order in the coupling constant, the density matri- 
ces of the system and of the environment factorize at all 
times. 

The behaviour of the effective decay rate, appearing in 
Eq. |J3}, defines the occurrence of the Zeno or anti-Zeno 



effect. We indicate with 7° the Markovian decay rate of 
the survival probability, as predicted by the Fermi Golden 
Rule. If a finite time r* such that 7^ (r*) = 7° exists, 
then for r < r* we have Jni^/ln < 1> i- e - the mea- 
surements hinder the decay (QZE). On the other hand, 
for r > t* we have 7^f (t)/t£ > 1 and the measurements 
enhance the decay (AZE) |l6j . 

For an initial Fock state \n) there exist two possible 
decay channels associated with the upward and down- 
ward transitions to the states |n -I- 1) and \n — 1), re- 
spectively. The probability that the oscillator leaves its 
initial state after a short time interval r can be written 
as P n (r) = P„ T (r) + P+(t), where PJ(t) and P^{t) are 
the probabilities that an upward or a downward tran- 
sition, respectively, has occurred in the interval of time 
< t < t. From Eq. (12) of Ref. 0], neglecting second 
order processes, one gets 



dt 



= -{(n+l)[A(t)-7(*)]+n[A(t)+7(*)]}frm(t),(4) 



where p nn {t) = (n\ps(t)\n). From Eq. (@J, noting that 
p nn (t) ~ Pnn(0), it is straightforward to see that the 
upward and downward transition probabilities in the in- 
terval < t < t are P^(t) = (n+1) £[A(t)-j(t)]dt and 
P 7 {(t) — n fy[A(t) +j(t)]dt, respectively. We note that 
the survival probability P n (T~) defined in Eq. © is given 
by P„(t) = 1 — Pn{i~). Assuming that r is small enough 
and keeping only the first two terms in the expansion of 
the exponential appearing in Eq. J3Jl one easily gets 



(2n + l) / A(t)dt 



j(t)dt 



(5) 



We notice that the quantity 7jf (t) can also be derived 
starting from the generalized master equation obtained 
in Ref. 0], applying the formalism to the case of the 
harmonic oscillator. By definition A(t) and "f(t), up to 
the second order in the coupling constant, are given by 

HI 

ft p 00 

A(<) = aJ^J(j)cotlXftw/2fc B P)cos(^i)cos(woii)dcj^i(6) 

f-t f-OO 

7(t) = a 2 J(<jj)sm(u;ti)sm(uJoti)dujdti, (7) 



with a microscopic dimensionless system-reservoir cou- 
pling constant. We note that, for high reservoir temper- 
atures T, A(t) > 7(4). Inserting Eqs. @-0 into Eq. © 
and carrying out the double time integration, one obtains 

7^)=r/jH {nfc-Jsincfa-T) + t»£ + Wc^w + t)} dw,(8) 



(±) 



with sinc(a;) = sin(x)/x, lo± — (u> ± loq)/2, and n 
a 2 [coth(hu/2k B T) (ra+l/2)±l]. In the limit t -> 00 one 
gets the Markovian value of the effective decay rate 7° = 
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(2n + 1)Am + 7m 5 with Am and 7m Markovian values 
of the diffusion and damping coefficients respectively. 

The quantity ruling the occurrence of either the QZE 
or the AZE is the ratio 



7l f (r) _ (2n + 1) ft A(t)dt - f Q 7 (t)di 



r[(2n + 1)A M - lM ] 



(9) 



It is worth underlining that, in general, this quantity 
depends on the initial state of the system, i.e. on the 
initial Fock state |n). For high T, however, since A(t) ^> 
7(t), Eq. (JSJ becomes independent of the initial state 



7^) „ / T A(t)dt 



7;; 



-A 



M 



(10) 



Equation IjlUI) approximates Eq. @ also for initial highly 
excited Fock states, i.e. n ^> 1. This equation, together 
with Eq. @ , establishes a connection between two funda- 
mental phenomena of quantum theory, namely the quan- 
tum Zeno effects and environment induced decoherence. 
This connection and its physical consequences, which will 
be carefully described in the rest of the paper, constitute 
our main result. We underline that, in deriving Eqs. 
(|lUp. no assumption on the form of the spectral density 
has been done. 

Equation l|10|) tells us that, for high T reservoirs and/or 
n>l the effective decay rate depends only on the diffu- 
sion coefficient A(i), which in turn describes EID. This 
yields a new physical explanation of the occurrence of ei- 
ther Zeno or anti-Zeno effects for QBM . The eigenstates 
of the quantum harmonic oscillator are highly nonclas- 
sical states. They are not localized, either in position 
or in momentum, and therefore they are very sensitive 
to environment induced decoherence. If the average EID 
in the interval between two measurements, quantified by 
J Q A(t)dt/r, is smaller than the Markovian one, quanti- 
fied by Am, then when the system "restarts" the evolu- 
tion after each measurement, the effect of EID is again 
less than in the Markovian case. Essentially the mea- 
surements force the system to experience repeatedly an 
effective EID which is less strong than the Markovian one. 
In this case the QZE occurs. On the contrary, whenever 
the average increase of EID in the interval between two 
measurements is greater than its Markovian value, then 
the measurements force the system to experience always 
a stronger decoherence, and hence the system decay is 
accelerated (AZE). 

An important consequence of Eq. illUf) is that the 
quantum-classical border can be manipulated by means 
of the QZE and of the AZE. When the QZE occurs, in- 
deed, an initial delocalized state of the harmonic oscilla- 
tor such as its energy eigenstate remains delocalized for 
longer times, hence the QZE "moves forward in time" 
the quantum-classical border. The opposite situation 
happens when the AZE occurs. 

The analytic expression of the diffusion coefficient al- 
lows us to single out the relevant system and reservoir 
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FIG. 1: (a) Ratio between the effective decay rate 7 ^f (r) and 
the Markovian decay rate 70 for high T. Solid line r = 0.5, 
dotted line r = 1, dashed line r = 10. (b) Ratio between 
the diffusion coefficient A(r) and its Markovian value Am for 
high T. Solid line r = 0.5, dotted line r = 1, dashed line 
r = 10. (c) Ratio between the effective decay rate 7 jf (t) and 
the Markovian decay rate 70 for T ~ and n S> 1. Solid line 
r = 0.1, dotted line r — 0.5, dashed line r = 10. (d) Ratio 
between the diffusion coefficient A(r) and its Markovian value 
Am for T ~ 0. Solid line r = 0.1, dotted line r = 0.5, dashed 
line r — 10. A comparison between (a) and (c) shows the 
temperature dependence of the Zeno-anti-Zeno crossover. 



parameters ruling the crossover between the QZE and 
the AZE. The time r* at which 7 Z (r*) = 70, when it ex- 
ists and is finite, can be seen as a transition time between 
Zeno and anti-Zeno phenomena |l6| . Our analysis shows 
that there exist two other relevant parameters, namely 
the ratio r = lu c /ujq, quantifying the asymmetry of the 
spectral distribution, and the ratio ksT/HuJo. For a fixed 
time interval r, a change in the values of r and T may 
lead to a passage from Zeno dynamics to anti-Zeno dy- 
namics and vice versa. 

For high T, it is easy to prove that a Zeno-anti-Zeno 
crossover exists only for values of r < 1, i.e. for luq > u> c , 
as it is shown in Fig^ (a). This can be understood in 
terms of EID by looking at the short time dynamics of 
A(t) in the time interval < t < r, i.e. before the first 
measurement is performed (see Fig. 0(b)). 

The situation changes drastically for the case of a zero- 
T reservoir, characterized by an asymmetric spectral den- 
sity. For n 3> 1, indeed, in contrast to the high T 
case, the Zeno-anti-Zeno crossover exists also for values 
of r S> 1. i.e. in the case luq <ti uj c . The region in which 
only Zeno dynamics may occur now appears at the edge 
of the spectral density function, i.e. for ujq ~ lu c . The 
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reason of such a different dynamics for small values of loq 
stems from the strong decoherence, showing up at low 
temperatures not only when r -C 1 (as in the high T 
case), but also for r ^> 1 (See Figs. Q ( c )-(d)). Summa- 
rizing, the occurrence of the Zeno or anti-Zeno effects is 
directly related to the absence or presence, respectively, 
of the so-called initial jolt of the diffusion coefficient A(t) 
0, which is the signal of strong initial decoherence. 

Another interesting aspect stemming from our re- 
sults concerns the Zcno-anti Zeno crossover for an initial 
ground state (n = 0) of the system oscillator. The high 
T behaviour is given by Eq. {TUJI [See Fig. [I] (a)- (b)]. For 
very low reservoir temperatures (T ~ 0), Am — 7m and 
therefore the denominator of Eq. (0 approaches zero, im- 
plying that J Z (t)/jq 3> 1 always, i.e. the measurements 
always enhance the decay (AZE). Summarizing, in this 
case, by changing the reservoir temperature, e.g. starting 
from high T and lowering the temperature, one observes 
a passage from the situation depicted in Fig. (a) to a 
situation in which only the AZE is practically observable. 

The possibility of controlling both the environment 
and the system-environment coupling would allow one to 
monitor the transition from Zeno to anti-Zeno dynam- 
ics. The use of artificial controllable engineered environ- 
ments has been recently demonstrated for single trapped 
ions In Ref. 19] it has been shown that the en- 

gineered amplitude reservoir realized by applying noisy 
electric fields to the trap electrodes in |l2( can be used 
to simulate quantum Brownian motion and to reveal the 
quadratic short time dynamics. Shuttering these noisy 
electric fields one could model a fast switch off-on of the 
environment. Actually, when the noise is off, the reser- 
voir simply does not exist anymore. The action of the 
sudden switch off-on of the environment may be seen as 
a physical implementation of the operation of trace over 
the reservoir degrees of freedom. The operation of trace 
is a typical example of a non-selective measurement (see, 
e.g., Ref. 0, p. 321). Hence a succession of short switch 
off-on periods, realized by shuttering the engineered ap- 
plied noise, would induce Zeno or anti-Zeno dynamics 
depending on the value of the system and reservoir pa- 
rameters and of the shuttering period. This is the core 
idea for monitoring the Zeno-anti-Zeno crossover with 
trapped ions. 

It is worth underlining that since the measurements 
causing the acceleration or inhibition of the decay, imple- 
mented by a fast switch off-on of the noisy electric field, 
are non-selective, they do not disturb the vibrational 
state of the ion. Therefore, by using the set up of Ref. 
[l^ . a comparison between the population P,[ N \t) of 
the initial vibrational state (e.g. the vibrational ground 
state \n = 0)) in presence of shuttered noise (with N 
the number of switching off-on periods), and the popula- 
tion P n (t) in presence of un- shuttered noise, would show 
that Pi N) (t) > P n (t) (QZE) or p[ N \t) < P n (t) (AZE) 



depending on the choice of the parameters. 

All the parameters ojq, ui c , T and r, driving the Zeno- 
anti-Zeno crossover, may be varied in the experiments. 
Although the value of ojq may be modified only within a 
certain range and under certain constrains, the modifica- 
tion of both lo c and T may be obtained by simply filtering 
the applied noise and varying the noise fluctuations, re- 
spectively |l9|. Since all the parameters ruling the Zeno 
- anti-Zeno crossover are easily adjustable in the experi- 
ments, its observation in the trapped ion context should 
be already in the grasp of the experimentalists. 
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